Light echoes reveal an unexpectedly cool 77 Carinae dur- 
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r] Carinae is one of the most massive binary stars in the Milky WajSOl It became the 
second-brightest star in the sky during its mid- 19th century "Great Eruption," but then 
faded from view (with only naked-eye estimates of brightnesses). Its eruption is unique 
among known astronomical transients in that it exceeded the Eddington luminosity limit 
for 10 years. Because it is only 2.3 kpc away, spatially resolved studies of the nebula have 
constrained the ejected mass and velocity, indicating that in its 19th century eruption, r] Car 
ejected more than 10 M in an event that had 10% of the energy of a typical core-collapse 
supernova 5 6 without destroying the star. Here we report the discovery of light echoes of 
r] Carinae which appear to be from the 1838-1858 Great Eruption. Spectra of these light 
echoes show only absorption lines, which are blueshifted by —210 km s -1 , in good agreement 
with predicted expansion speeds 6 . The light-echo spectra correlate best with those of G2-G5 
supergiant spectra, which have effective temperatures of ~5000 K. In contrast to the class of 
extragalactic outbursts assumed to be analogs of eta Carinae 's Great EruptiorP^, the effec- 
tive temperature of its outburst is significantly cooler than allowed by standard opaque wind 
models^. This indicates that other physical mechanisms like an energetic blast wave may 
have triggered and influenced the eruption. 
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rj Car-like giant eruptions of Luminous Blue Variables (LBVs) are characterized by signifi- 
cant mass-loss and an increase in luminosity by several magnitudes 8 9 . Traditionally, it was thought 
that this increase in luminosity drives a dense wind, producing an optically-thick cooler pseudo- 
photosphere with a minimum effective temperature of 7000 K and an F-type spectrum^. Within 
this model, 77 Car has been considered the prototype of the so-called "supernova imposters"^^. 

We obtained images in proximity to 77 Car (Figure [T} that, when differenced, show a rich set 
of light echoes. The largest interval between our images is 8 years. We have also found similar 
echo candidates at other positions, which we are currently monitoring. The large brightening and 
long duration point to the Great Eruption as the source of the light echoes. We have also obtained 
a composite light curve in the SDSS i filter of the light echoes (see Fig. [2]), showing a slow decline 
of several tenths of a magnitude over half a year. This light curve is most consistent with the 
historical observationPof a peak in 1843, part of the 1838-1858 Great Eruption, although further 
observations are necessary to be certain (see the Supplementary Information, SI). 

Spectra of the light echoes (see Figure [3]) show only absorption lines characteristic of cool 
stellar photospheres, but no evidence of emission lines. In particular the Ca II infrared (IR) triplet 
is only in absorption. Because of bright ambient nebular emission, it is difficult to determine if 
there is any Ha emission from 77 Car itself, but in any case it must be weak if present. By cross- 
correlating each of our 77 Car echo spectra with the UVES spectral library 14 (see Figure S2 and S3) 
we find best agreement with supergiant spectral types in the range of G2-G5, with an effective 
temperature of ~5000 K. Spectral types of F7 or earlier are ruled out by our analysis (see SI for 
more details). 

Doppler shifts of absorption features in the echo spectra provide direct information about the 
outflow speeds during the eruption. The Ca II IR triplet absorption features in the spectrum are 
noticeably blueshifted (see Figure S2). By cross-correlation with G-typeP^ templates, we determine 
velocities of —202 ± 9, —210 ± 14, and —237 ± 17 km s _1 for our three individual spectra and an 
average velocity of —210 ± 30 km s -1 , which includes an uncertainty for the dust sheet velocity. 

The bipolar nature of the Homunculus shows that the 77 Car Great Eruption was strongly 
aspherical. It was previously predicted that the outflow speeds one would derive from spectra of 
77 Car in outburst, looking at the poles and equator of the double lobes, would be ~ —650 km s 1 
and —40-100 km s _1 , respectively 16 (outflow speeds near the equator have a steep latitude depen- 
dence). The light echo we investigate in this paper arises from latitudes near the equator of r/ Car 
(see Figure SI), and the measured blue-shifted velocity of —210 ± 30 km s _1 is in good agreement 
with expansion speeds within ±20° of the equatorial plane. We also find a strong asymmetry in 
the Ca II IR triplet, extending to a velocity of —850 km s _1 - well below the speed of the fastest 
polar ejecta found previously 6 , but is in good agreement with speeds observed in the blast wave 
at lower latitudes 6 . Future observations of light echoes viewing the 77 Car eruption from different 
directions, in particular from the poles, has the potential to observe these very high-velocity ejecta 
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and other asymmetries. 

A characteristic of LBV outbursts is their transition from a hot quiescent state to a cooler 
outburst state, although this feature is less well observed for the giant eruptions (see Figure [4]). 
Two potential models for LBV outbursts involve either an opaque stellar wind driven by an in- 
crease in luminosity, or a hydrodynamic explosion. The traditional mechanism for 77 Car-like 
giant eruptions has been that an unexplained increase in luminosity drives a denser wind, so that 
an optically thick pseudo-photosphere forms at a layer much larger and cooler than the hydrostatic 
stellar surface 13 . This model predicts a minimum effective temperature of 7000 K, resembling A or 
F-type supergiants^EZl, since the wind opacity depends on the temperature (see Figure|4]). They ev- 
idently occur as a massive star attempts to evolve redward and encounter the Humphreys-Davidson 
Limit 13 , beyond which no stable stars are observed. 

Surprisingly, our G-type light-echo spectrum of the 77 Car Great Eruption is inconsistent 
with expectations of an opaque-wind model 13 (see Figure [i}. This model also fails to explain the 
high 10 50 erg kinetic energy 5 and the presence of a fast blast wave at large radii 6 . Instead, these 
observations point toward a hydrodynamic explosion at least influencing the Great Eruption^EIl. 

The first visual spectroscopic observations of rj Car around 1870 showed emission lines 18 19 . 
A photographic spectrogram obtained during its Lesser Eruption circa 189C 20 21 resembles an F- 
type supergiant blueshifted by —200 km s _1 , with moderate hydrogen P Cygni profiles, which is as 
expected in the opaque- wind modeP^ The difference between the 1890 and our light-echo spectra 
of the Great Eruption is therefore quite striking, indicating that two distinct physical processes may 
have been involved for two outbursts of the same object. However, the 1890 event also produced a 
mass ejection, the Little Homunculus, with the same axial symmetry as the Great Eruption 22 , albeit 
of a much smaller amount. 

LBV giant eruptions are rare, and only have been recorded in our Galaxy in the last 400 
years, the Great Eruption of 77 Car and P Cygni's giant eruption in the 17th century. Because 
of their considerable intrinsic brightness just below the luminosity of faint core-collapse SNe, 
about two dozen giant eruption candidates, so-called SN imposters since they have been often 
mistaken for SNe, have been found in various extragalactic SN searches 7 12 . Typically, the hotter 
SN imposters have steep blue continua, stronger and broader Balmer lines, and relatively weak 
absorption, whereas the cooler ones tend to have redder continua, weaker and narrower Balmer 
lines, strong [Ca II] and Ca II emission, deeper P Cygni absorption features, and in some cases 
stronger absorption spectra similar to F-type supergiants 23 . However, the r\ Car Great-Eruption 
light-echo spectrum is quite different. Its spectral type is G2-G5, significantly later than all other 
SN imposters at peak. Furthermore, the Ca II IR triplet lines are only in absorption. For the 
extreme mass-loss rates required in 77 Car's Great Eruption, another process must give rise to the 
apparent temperature. 
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rj Car's Great Eruption has been considered the prototype of the extragalactic SN imposters or 
7] Car analogues, even though it is actually an extreme case in terms of radiated energy (10 49,3 erg), 
kinetic energy (>10 50 erg), and its decade-long duration^. The spectra of the light echo indicates 
now that it is not only extreme, but a different and rather unique object. It is difficult to see 
how strong emission lines could be avoided in an opaque wind where the continuum photosphere 
is determined by a change in opacity, and its temperature and broad absorption lines are more 
consistent with the opaque cooling photosphere of an explosion. The cause that triggered such 
an explosion and the mass-loss without destroying the star is still unknown, but predictions from 
future radiative transfer simulations trying to explain t] Car and its Great Eruption can now be 
matched to these spectral observations. Other alternative models that were proposed, e.g. the 
ones that use mass accretion from the companion star during periapsis passage as a trigger for the 
eruption 12 , can be either verified or dismissed. I-H-I 
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SI Supplementary Information 



Three dimensional orientation The scattered-light path is shown in Figure [ST] The 3D-plot 



shows that with this light echo we see the eruption of rj Car from a viewing angle perpendicular to 
the principal axis of the Homunculus Nebula. It will be very interesting to compare spectra seen 
from this viewing angle with spectra from light echoes that view the eruption from the poles of the 
bipolar outflow. 



Evidence that these are light echoes of the mid- 19th century Great Eruption of rj Car The 

angular separation between the center of our light-echo images and rj Car is ~0.5 degrees, as shown 
in the left panel of Figure 1, in which our pointing is indicated with a white box. In the middle 
panel of this figure, our three epochs of SDSS /-band imaging obtained with the 8k x 8k Mosaic 
imager on the CTIO 4-m Blanco telescope are displayed. The images, each with an exposure 
of 160 s, were obtained on 2003 Mar 10 (A), 2010 May 10 (5), and 2011 Feb 6 (C). Difference 
images, C — A and C — B, are also shown (upper and middle right panels, respectively). In them, 
we notice excess flux from the first (black) and second epochs (white), which we interpret as light 
echoes from rj Car. Sample light-echo positions are indicated with blue (epoch A) and red (epoch 
B and C) arrows. 



Comparing the first- and second-epoch images, we can determine the direction of the appar- 
ent motion of the light echo and infer the direction to the illuminating source. Most of the excess 
flux in the first-epoch image is toward the northern part of the image (blue arrows in Figure 1, 
whereas most of the excess flux in the second-epoch image is in the southern portion (red arrows 
in Figure 1 . If the light echoes arise from r\ Car, this is exactly the directional sense expected. We 
have also found similar regions of excess flux at other positions in proximity to rj Car, with ap- 
parent motions consistent with rj Car being the source of the outburst light, which we will discuss 
in a follow-up paper. The vector method we have used previously to determine the origin of the 
light echoes for SNR 0509-67.5, Cas A, and Tycho^ED was applied here and, not surprisingly, 
we conclude that these excess fluxes are most likely light echoes from a dramatic brightening of 
rj Car. 



The flux profile of a light echo, which is the cut through the light echo along the axis pointing 
toward the source event, is the projected light curve, stretched or compressed depending on the 
inclination of the scattering dust filament, and convolved with the effects of the dust width and 
the seeing 26 . In pathological cases where the inclination of the dust filament is very unfavorable 
or where the dust filament is significantly thicker than any yet observed, a useful correspondence 
may be losP^. For light echoes from sources in our Galaxy scattered by dust with typical widths 
and inclinations, we expect that the apparent motion of the light echo will move at a rate of one 
projected echo width between two epochs separated by twice the outburst event duration. For 
known galactic SNe light echoes the timescale of the brightest phase is of the order of a few 
months, so that an imaging separation time of 1 year is sufficient for there to be no overlap (on 
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the sky) of light-echo featuresS3EH. For a several-decade duration outburst, like that of 77 Car, the 
difference image C — B has only a time difference of 1 year, and can only reveal the small boundary 
fractions of the light echo in C, whose more complete extent can be seen in C — A which has a time 
difference of 8 years. We conclude that the duration of the event illuminating these light echoes 
is significantly longer than one year, consistent with our conclusion that these are echoes of the 
20-year long Great Eruption of 77 Car. 



Determination of spectral type We compare the light-echo spectra to the compilation of super- 
giant spectra in the UVES atlas^l (see Table S2). The upper left panel of Figure S2 shows the 
light-echo spectra, flattened by dividing by low-order continuum fits, and the UVES supergiant 
sequence. The instrument configuration for the EC1A spectrum was designed to optimize the S/N 
ratio in the wavelength region of the Ca II IR triplet, resulting in a lower S/N ratio for A < 6000 A. 
For this reason the EC1A spectrum is absent from the upper left panel of Figure|S2} The light-echo 
spectra correlate very well with late-F and G-type stars, in particular the Mg b lines and line blends 
at 5270 A and 6497 A. Both earlier- and later-type spectra show significantly fewer similarities. 



We calculate the cross-correlation parameter r between the light echo and the UVES spectra 
using the IRAF routine xcsao in the wavelength range 5050-6500 A (see left panel of Figure S3 ). 
We exclude wavelength ranges contaminated by fore/background emission lines (see Table S4). 
We find that extending the wavelength region to redder wavelengths (e.g., 5050-7500 A) decreases 
the correlation, and in particular the correlation differences between the different spectral types. 
The reasons for a decrease in correlation are threefold: (1) the spectral differences between F, G, 
and K stars are stronger in the bluer wavelength ranges, and therefore including redder wavelengths 
dilutes the discrimination power, (2) strong fore/background emission lines like H a , N II, and 
S II introduce large discontinuities in the wavelength coverage, and (3) the chip gaps fall into the 
6500-7500 A wavelength range, further adding to the discontinuities. For the above reasons we 
determine the 5050-6500 A wavelength range to be optimal for the correlation analysis. 



To quantitatively estimate the temperatures of the supergiants that correlate best, we smooth 
r(T e ff) with a Gaussian of width a = 300 K (see lines in left panel of Figure S3 ). We find that the 
temperature with maximum correlation is 5210 K and 4950K for EC1B and EC1C, respectively. 
The right panel of Figure [S3] shows the probability density functions (PDFs) of the temperature, 
T e ff, of the supergiant spectral template having the best smoothed correlation, r(T e ff), with each 
light-echo spectrum. The PDFs are computed by bootstrap resampling of the distribution 10 5 
times. The 95% confidence intervals of the best-matching temperatures are 4850-5550 K (spectral 
type G0-G5) and 4450-5400 K (spectral type G0-K1) for EC1B and EC1C, respectively. Using 
an extended wavelength range of 5050-7500 A, the 95% confidence intervals are 3700-6300 K 
(F8-K5) and 3450-5100 K (G2-M1) for EC1B and EC1C, respectively. For the reasons mentioned 
above, the extended wavelength range results in a spectral type that is less constrained compared to 
the 5050-6500 A range. We conclude that the light-echo spectra agree best with supergiant spectral 
types in the range of G2-G5, with an effective temperature of ~5000 K. We conservatively exclude 
spectral types of F7 and earlier. 
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Unfortunately, the UVES spectra have a gap in coverage around the Ca II IR triplet, and 
therefore we use a compilation of supergiant spectra 15 covering the Ca II triplet wavelength range 
(see Table S3). The lower- left panel of Figure S2 compares the observed light-echo spectrum 
with this atlas. For early-type stars up to GO, the wavelength range redward of 8300 A is mainly 
dominated by the H Paschen series, which contaminates the Ca II triplet but gradually disappears 
as the effective temperature decreases 15 . Note that the observed light-echo spectra do not show any 
clear signs of the H Paschen lines; in particular, the lines at 8409 A and 8498 A are at most very 
weak. The earliest spectral types in agreement with such weak or non-existent H-Paschen lines are 
late F-types, in excellent agreement with the spectral correlation results in the wavelength range 
5050-6500 A. 



Visual observers noted a reddish or "ruddy" color during r\ Car's Great Eruption, e.g. de- 
scribed by Herschel as "redder than Arcturus'^. These observations point to a temperature of 
T e ff < 4500 K, lower than the temperature we infer from the spectral lines. It is plausible that 
reddening by grains forming in the eruption cause the apparent color. We note, however, that the 
temperature inferred from the spectral type is much more reliable than an apparent color, in partic- 
ular one done visually, since it is not influenced by unknown amounts of reddening caused by new 
dust formation along the line of sight. 



The lower-right panel of Figure [S2] shows velocities determined from cross-correlation from 
spectra of different spectral types. The contribution to the derived velocity caused by the motion 
of the reflecting dust sheet, i.e. the moving-mirror effect, is likely less than 30 km s -1 given 
the relatively low dominant expansion speeds of cool gas in the Carina Nebula H II region 29 . This 
systematic error could be constrained better with more spectra of light echoes located close enough 
to each other that they probe similar viewing angles, but different enough that the scattering dust 
is parsecs apart and thus has independent velocities. 



The most similar SN impostor The closest SN impostor in temperature is UGC 2773 OT2009-1, 
which is also dominated by a forest of absorption lines, similar to an F-type supergianPZMl. How- 
ever, the Ca II IR triplet has a P Cygni profile with a strong emission component 28 . Interestingly, 
UGC 2773 OT2009-1 is also one of the few examples where the outburst has persisted for years 
and in fact still continued 23 ! 27 ! 28 !. 



Ha line The interpretation of the Ha emission line is difficult, since the fore/background emission 
line subtraction is imperfect. The left panel of Figure 3 shows the Ha line of the three observed 
light-echo spectra. Both EC1B and EC1C show narrow emission at zero velocity, most likely 
caused by the incompletely subtracted fore/background emission. However, they also show asym- 
metry at the wavelength for which we would expect Ha emission if the hydrogen has the same 
blueshifted velocity of -210 km s _1 , indicated with the red line. If this is indeed a blueshifted Ha 
from the eruption, it is a narrow one unlike many of the SN impostors, but similar to the coolest 
ones like UGC 2773 OT2009- lEDH . 
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Epoch of Spectra The slit position of the light echo spectrum is shown in the bottom right panel 
of Figure 1 . Because the 77 Car Great Eruption lasted over a decade, it is not straightforward to 
determine the epochs probed by the light-echo spectra. In order to attempt this, we compare the 
historical light curve 4 (see Figure 2) with time variation of the light echoes. 

We summarize here an accounP of the photometric history of the Great Eruption of 77 Car. 
John Herschel recorded a brightening of ~1 mag in less than two weeks near the end of 1837. 
Over the next few months, 77 Car slowly faded again. Unfortunately, there is a gap in brightness 
estimates for the period between late 1838 and 1841 . There was another brightening in early 1843, 
and then another episode where 77 Car again faded back to typical non-outburst brightness. Late in 

1844, i] Car's greatest recorded brightening episode began, with peak brightness occurring in early 

1845. Thereafter a slow decline began which lasted for the next decade. 

We can compare 77 Car's historical light curve to the light curve we derived from the light 
echoes: at a given epoch and position on the sky, the light-echo flux is the flux of the source event 
integrated over a range of epochs. The range of epochs represented in the instantaneous flux at a 
given pointing depends on both astrophysical (dust width and inclination) and observational (see- 
ing) factors, and can be represented by a near-Gaussian window functional For the light-echo 
systems we have characterized around the Cas A and Tycho remnants, the full width at half maxi- 
mum was on the order of 25-120 days 26 . Thus if the flux at the same spatial position is measured 
at different epochs, the brightness recorded will be the true outburst light curve convolved with the 
window function. 

With four more epochs of the light echo obtained at the Faulkes Telescope South (FTS) 2-m 
telescope, we generated light curves of 77 Car at 17 positions along the spectroscopic slit shown 
in Figure 1. Since we have images from 2003 to mid-2011, these light curves span ~8 years. 
After an overall normalization to each of these light curves to match each other, the differences 
in light-curve profiles are remarkably small, with a typical standard deviation in flux of 2-4% (see 
Table S5). We therefore average these light curves into a single light curve. We show this light 
curve in Figure 2, shifted by 174.2 yrs (green circles), 167.95 yrs (red circles), and 166.28 yrs 
(blue circles) in order to match the 1838, 1843, and 1845 outbursts, respectively. The earliest 
epoch (2003) needs to be considered an upper limit, since that epoch was used as the template 
for the difference imaging, and we determined the light-echo flux at that epoch as the difference 
between the flux in the template image and the flux at a different, apparently empty position in the 
same image. In other words, there is an undetermined flux zeropoint when an on-frame reference 
location may contain light echo flux. It is also important to note that the light echo light curve is 
determined from images in SDSS i, which can introduces systematics when comparing it to the 
visual magnitude of the historical light curve, in particular since LBV's change their colors during 
outbursts. 

The light-echo light curve we derive from our images clearly reveals brightening of 2 mag- 
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nitudes or more within an interval of 8 years. The Lesser Eruption of rj Car from 1887 to 1896 
caused r\ Car to brighten by roughly 1 magnitude only. Unless the Lesser Eruption was extremely 
asymmetric such that it increased by an additional 1 magnitude or more in the direction of the light 
echo, the Lesser Eruption is excluded as the source of the light echoes. We also note that most of 
the brightening occurs between the second and third epochs in a 9 month time- span, indicating that 
the peak was at some time between these epochs. 

We conclude that there are three remaining eruption epochs to which our observed light curve 
could be assigned: 1) the 1838 peak, 2) the 1843 peak, and 3) the 1845 peak. A fourth possibility 
is an even earlier, unrecorded brightening. Our ability to decide among the different epochs relies 
on the time differences between our images and spectra. The second- and third-epoch pair have 
only 9 months between them and therefore provide strong discrimination for brightenings. The 
spectroscopic observations were done shortly after the third-epoch images, corresponding to <9 
months after a brightness maximum. 

We now discuss the individual epochs and the likelihood that they correspond to the bright- 
ening observed in our images: 

• 1838 peak: The first two epochs agree very well with the historical light curve. The histori- 
cal light curve seems to decline faster after the peak at the beginning of 1838. Unfortunately, 
the historical data that constrain this decline are few and have large uncertainties. In addition, 
if the width of the scattering dust filament is very large, then the window function can span 
a wide range of epochs, which would cause the light-echo light curve to have a shallower 
decline then the event light curve. 

• 1843 peak: The first two epochs again agree very well. The decline in brightness is in very 
good agreement with what would be expected from the historical light curve. 

• 1845 peak: The first two epochs also agree very well with the 1 845 peak. However, the cur- 
rent decline in brightness mag yr _1 ) is considerably faster than the typical decline of 
~0.1mag yr -1 observed after 1845. The light-echo light curve can only be put in agreement 
with the historical light curve if there was a fast, unobserved decline from the peak in 1845. 

We conclude that the most likely source of the light echoes for which we have spectra is 
the 1843 outburst, shortly after its peak. Assuming typical dust properties as observed for light 
echoes of other Galactic sources 26 \ the spectra is the light-curve-weighted average over an epoch 
range of 1-4 months. In the future, we will be able to further constrain the epoch range with 
additional imaging data. The two peaks in 1838 and 1845 cannot yet be excluded completely, but 
continued observations of the light echoes at monthly intervals should discriminate the possibilities 
within a year. This work, combined with targeted new observations and other contemporary images 
which may become known, will likely provide a detailed photometric time-series of the Great 
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Eruption with the probability of a useful companion spectral time-series. The recovery of valuable 
astrophysical observational data from the pre-imaging era is another powerful illustration of the 
power of light-echo observations 
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Figure 1: rj Car light echoes. The left panel shows the positions of rj Car and our images (white 
box), plotted on an image in the light of 3 different emission lines: oxygen (blue), hydrogen 
(green), and sulfur (red) (credit: Nathan Smith, University of Arizona/NOAO/AURA/NSF). The 
middle panels show the images obtained with the CTIO 4-m Blanco telescope of a region ~0.5 de- 
grees to the south of rj Car at epochs 2003 March 10 (A), 2010 May 10 (B), and 201 1 Feb 6 (C), 
from top to bottom. The right panels show the difference images C-A and C-B at the top and mid- 
dle, respectively. Example light-echo positions are indicated with blue (epoch A) and red (epoch B 
and C) arrows. The bottom right panel shows a zoom of the spectrograph slit, indicated with a blue 
line. For all panels north is up and east is to the left. Applying the vector method that previously al- 
lowed us to identify the source of the light echoes from the supernovae (SNe) that produced the SN 
remnants 0509—67.5, Cas A, and Tychd^O^, we find that a dramatic brightening of r\ Car must be 
the origin. In these echoes, unlike those of Galactic SNe 26 , there is still significant spatial overlap 
even at separations of one light-year, suggesting that the duration of the event causing them must 
be significantly longer than one year. We also see brightening of 2 magnitudes or more within 8 
years. Thus, the Lesser Eruption from 1887 to 1896, which brightened by only a magnitude, is 
excluded as the source. 
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Figure 2: Historical and light-echo lightcurve of r\ Car. The historical light curveP in visual ap- 
parent magnitudes is shown with black circles and lines, with error bars indicating approximate 
uncertainties in these eye estimates. Light echo brightnesses (SDSS i; error bars are the s.d.) from 
our eight modern images spanning ~8 years are displayed shifted by 174.2 yrs (green circles), 
167.95 yrs (red circles), and 166.28 yrs (blue circles), in order to illustrate the best-matching time 
delays for the 1838, 1843, and 1845 outbursts, respectively. The first epoch is an upper limit indi- 
cated with an arrow. The upper panel shows the full time range of the Great Eruption and therefore 
shows all three potential matches, whereas the lower panels show the brightnesses from seven of 
our eight modern epochs in a magnified time period around each peak. 
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Figure 3: Light echo spectra of rj Car's Great Eruption. The three optical low-resolution 
spectra of the light echo (black lines) were taken at J2000 position RA=10:44: 12.127 and 
Decl.=-60: 16:01. 69 in March and April 2011 obtained at the Magellan I 6.5-m and du Pont 2.5- 
m telescopes of the Las Campanas Observatory, Chile. A log of the spectroscopic observations 
and details of the spectra is presented in Table SI. The slit positions differ only slightly in slit 
angle. The spectra were reduced using standard techniques and then wavelength-calibrated using 
observations of an HeNeAr lamp. The wavelength calibration was checked and corrected using 
night-sky emission lines, especially O I 5577A, and OH lines in the red part of the spectrum. We 
flux-calibrated the spectra using a flux standard observed the same night as the science observa- 
tions. The left panel shows the spectra from 5000 to 9000A. The spectra are not corrected for 
reddening nor for the blue-ward scattering by the dust. The blue lines show for comparison spectra 
of three examples of SN imposters: SN 1997bs, SN 2009ip, and UGC 2773-OT1. The right panel 
shows the Ha and [N II] emission lines. Note that the background emission-line subtraction is 
incomplete since they spatially vary. Also, for EC1A Ha is at the edge of the chip and is therefore 
uncertain. 
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Figure 4: HR diagram with LBVs and 77 Car. AdaptatiorP-^ of an HR diagram showing LBVs, 
related hypergiant stars, and the peak luminosities of LBV-like eruptions. The grey bands denote 
the typical locations of LBVs in quiescence and in S Doradus excursions. Temperatures for the 
Great Eruption and 1 890 eruption of 77 Car are based on the echo spectra presented here and the 
F-type spectrum of the 1890 evenP^ respectively. The temperature of 10,000 K for SN 2009ip 
is based on the observed continuum shape, but this is only a lower limit because of the possible 
effects of circumstellar or host galaxy reddening 27 . Because of the presence of He I lines in the 
spectrum, the true temperature is probably much hotter. The 8500 K temperature of UGC2773-OT 
is indicated by the F-type absorption features in its spectrum, and this temperature is relatively 
independent of reddening! 27 ! 28 !. 
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Figure SI: 3D-plot of the light path. North is toward the positive y-axis (up), east is toward the 
negative x-axis (left), and the positive z-axis points toward the observer with the origin at 77 Car. 
The red and brown circles indicate 77 Car and scattering dust, respectively. The parabolic relation 
between the spatial parameters of the scattering dust and the time since outburst is described by 
the well-known light echo equation® Assuming a time since outburst of 169 years, and a distance 
of 7660 light-years^, we find that the scattering dust is at a position (x,y,z)=(l4,-7&,-66) in light- 
years. The black lines show the path of the light scattering from the light-echo-producing dust 
concentrations. The top-right panel shows an HST image (credit: Nathan Smith / Jon Morse / 
NASA) of 77 Car, which is surrounded by expanding lobes of gas denoted as the Homunculus 
Nebula visible to the NW and SE. These lobes were created by the 1838-1858 Great Eruption. We 
indicate the lobes in the 3D plots with the two cones (Note that the size of the cones is not to scale). 
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Figure S2: Spectral type comparison. The stellar comparison spectra are convolved with a Gaus- 
sian of FWHM 7 A so that their resolution matches those of EC1B and EC1C. Similarly, we 
convolve the spectrum EC1A, which has an original resolution of 4A, with a Gaussian of FWHM 
5.7 A. We then divide the spectra by the low-order continuum which we determine by convolving 
the spectra with a Gaussian of FWHM 200 A. The upper left panel compares two of the observed 
light echo spectra to a selection of UVES supergiants 14 in the 5060-5500 A wavelength range (We 
do not show spectrum EC1A since its S/N ratio in this blue wavelength range is too low). The light 
echo spectra correlate very well with late-F and G-type stars, in particular the Mg b lines, and the 
Ca I, Fe I, Ti I, Cr I blend at 5270 A. The upper right panel shows the cross-correlation between 
the light echo and the UVES spectra in the wavelength range of 5050-6500 A as calculated by the 
IRAF routine xcsao. Since the UVES spectra have a gap around the Ca II IR triplet, we use the 
Ca II triplet spectral library^ in that wavelength range, as shown in the bottom left panel. The 
bottom right panel shows the blue-shifted velocity determined by xcsao from the Ca II triplet 
with respect to the effective temperature^ of the supergiant template spectra. 
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Figure S3: Effective temperature of the light echo spectrum. The left panel shows the cross- 
correlation between the light echo and the UVES spectra as calculated by the IRAF routine xcsao 
for EC IB (red squares) and EC1C (blue diamonds). The lines indicate r(T) smoothed with a 
Gaussian of width a = 300 K, peaking at 5210 K and 4950 K for EC1B (red) and EC1C (blue), 
respectively. The right panel shows the PDF of the best-correlating supergiant temperature, de- 
termined by bootstrap resampling the r distribution 10 5 times. The 95% temperature confidence 
intervals are 4850-5550 K and 4450-5400 K for EC1B and EC1C, respectively. 
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Star Name 


HD 




Spectral Type 




N Car 


HD 


47306 


AO II 


9500 


HR4541 


HD 


102878 


A2 lab 


9100 




HD 


34295 


A4 II 


8800 


n Vel 


HD 


74272 


A5 II 


8500 


y Car 


HD 


97534 


A6 la 


8400 


HR 3739 


HD 


81471 


A7 lab 


8300 


i Car 


HD 


80404 


A8 lb 


8200 




HD 


104111 


A8 II 


8200 


V399 Car 


HD 


90772 


A9 la 


8000 


HR 3496 


HD 


75276 


F2 lab 


7480 


b Vel 


HD 


74180 


F3 la 


7320 


HR 5024 


HD 


115778 


F4 II 


7160 


pPup 


HD 


67523 


F6 II 


6600 


BG Cru 


HD 


108968 


F7 Ib/ll 


6400 


HR 8470 


HD 


210848 


F7 II 


6400 


5CMa 


HD 


54605 


F8 lab 


6200 


V810Cen 


HD 


101947 


F9 lab 


5900 


7 1 Nor 


HD 


146143 


F9 la 


5900 


BB Sgr 


HD 


174383 


GO lb 


5500 


ERCar 


HD 


97082 


G1 lab/lb 


5300 




HD 


136537 


G2 II 


5100 


Crv 


HD 


109379 


G5 II 


4830 




HD 


125809 


G5/G6 lb 


4790 


HR 3673 


HD 


79698 


G6 II 


4750 


r Leo 


HD 


99648 


G8 lab 


4590 


d Cen 


HD 


1 1 7440 


G9 lb 


4500 


HR 3583 


HD 


77020 


G9 II 


4500 


HR 554 


HD 


11643 


K1 II 


4400 


e Peg 


HD 


206778 


K2 lb 


4300 


3 Cet 


HD 


225212 


K3 lab 


4000 


HR611 


HD 


12642 


K5 lab 


3750 


HR 2508 


HD 


49331 


M1 lab 


3450 


V528 Car 


HD 


95950 


M2 lb 


3350 


CRCir 


HD 


131217 


M2/M3 II 


3300 



Table S2: Supplementary Information: LIVES Supergiant Stellar Spectra Library 14 . We 
convert the spectral type into the effective temperature T eff using the spectral-type to 
temperature relation for supergiants ^ 
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Star Name 


HD 


Spectral Type 


To 




5 Per 


HD 13267 


B5 la 


13800 


13700 


v Cep 


HD 207260 


A2 la 


9100 


9100 


4> Cas 


HD 7927 


F0 la 


7425 


7800 


v Aql 


HD 182835 


F2 lb 


7350 


7480 


44 Cyg 


HD 195593 


F5 lab 


6600 


7000 


35 Cyg 


HD 193370 


F6 lb 


6200 


6600 


V440 Per 


HD 14662 


F7 lb 


5900 


6400 


HR 7008 


HD 172365 


F8 Ib-ll 


5500 


6200 




HD 18391 


GO la 


5500 


5500 


HR 7456 


HD 185018 


GO lb 


5550 


5500 


14 Per 


HD 16901 


GO Ib-ll 


5478 


5500 


HR 3459 


HD 74395 


G2 lab 


5250 


5100 


/3 Dra 


HD 159181 


G2 lab 


5250 


5100 




HD 187299 


G5 la 


5010 


4830 


/3Lep 


HD 36079 


G5 II 


5170 


4830 


£Pup 


HD 63700 


G6 la 


4990 


4750 


e Cet 


HD 13611 


G8 lab 


5040 


4590 




HD 12014 


KO lb 


5173 


4450 


CCep 


HD 210745 


K1.5 lb 


4500 


4350 


7 1 And 


HD 12533 


K3 Mb 


4383 


4000 


41 Gem 


HD 52005 


K4 lab 


4116 


3900 


V809 Cas 


HD 219978 


K4.5 lb 


4250 


3825 


a Ori 


HD 39801 


M2 lab 


3614 


3350 



Table S3: Supplementary Information: Ca II IR triplet Supergiant Stellar Spectra Library^. 
The temperature T is the effective temperature quoted by Cenarro et alM For consis- 
tency we also calculate T cS using the same stellar type-temperature relation 31 we used 
for the LIVES spectral library. 
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Line 




^max 


I 5577.4 


5572 


5582 


Na I D 5890.0, 5895.9 


5860 


5921 


I 6300.2 


6295 


6305 


I 6363.9 


6359 


6369 



Table S4: List of the fore/background emission lines which are excluded from the 5050- 
6500 Awavelength range of the correlation analysis. For each spectroscopic line, \ min 
and \ max define the initial and final wavelengths in A of the contaminating line region we 
exclude. 
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